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Abstract. Excessive fructose consumption was shown to have deleterious effects on the 
cardiovascular system, particularly as the metabolic syndrome. However, the degree by which 
alteration of each pathophysiological factor contributes to the morbidity associated with fructose 
consumption is not yet clear. We have developed a mathematical model to integrate and 
uniformly quantify pathophysiological features of the metabolic syndrome on a high fructose-fed 
dog model. A novel comprehensive measure for the syndrome severity (the “patholog”’) and a 
more intuitive measure of the insulin resistance are introduced. Alteration of hemodynamics, 
echocardiography and blood chemistry were determined in adult male mongrel dogs fed with 
60% isocaloric fructose or normal chow for 7 weeks. The diverse experimental data were 
transformed into comparable scores and a pre-Hilbert space of states constructed. In such a 
space one can quantify the severity of any combination of pathophysiological and genomic 
features and determine the global recovery resulting from a treatment. The model indicates 
increase of insulin resistance (new index proposed), systolic blood pressure, low-to-high density 
lipids ratio and angiotensin II as the major contributors to the excessive fructose morbidity. Our 
model provides the simplest, yet the most intuitive and comprehensive way to integrate data of a 
wide diversity in visualizing and quantifying a cardiovascular disease. 
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1. Introduction 

One major cause of obesity, type 2 diabetes, fatty liver and cardiovascular disease of 
millions in the Western world countries is the excessive sugar consumption, particularly 
in the form of high-fructose corn syrup [1-4]. Such consequences of the excessive 
fructose consumption are heralded by the metabolic syndrome [5-9], a combination of 
medical disorders reaching a prevalence of 35-39% in the United States [10]. According 
to the American Heart Association, the metabolic syndrome is characterized by: 1) 
increased waist circumference (men > 102 cm, women >90 cm), 2) elevated triglycerides 
(=150 mg/dL), 3) reduced high density lipids (men < 40 mg/dL, women (< 50 mg/dL), 4) 
elevated blood pressure (> 130/85 mmHg), 5) elevated fasting glucose (=> 100 mg/dL)I. 
High fructose diet results in hypertension [11], atherosclerosis [12], hyperurecimia [13], 
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endothelial dysregulation [14], dyslipidemia [15,16] consistent with the metabolic 
syndrome. 

Fructose is a monosaccharide naturally found in many fruits and vegetables. A 
major dietary source comes from the high-fructose corn syrup commonly added 
processed foods such as soft drinks, condiments, applesauce and even baby food [17]. In 
contrast to glucose which is metabolized in many organs, fructose is solely metabolized 
in the liver. Fructose metabolism bypasses phosphofructokinase in glycolysis and is not 
affected by insulin. Fructose actually induces insulin resistance [18-20], leading to a 
higher caloric intake by reducing the level of circulating leptin [21-23]. 


2. Materials and methods 

This report reanalyzes from a novel perspective the experimental data on fructose- 
fed dogs included in the PhD thesis of Dr. Kertowidjojo. Experimental protocols were 
approved by the Institutional Animal Care and Use Committee and respected the 
guidelines of the National Institutes of Health and the American Physiological Society. 

2.1. Animals 

Six adult mongrel male dogs were fed for seven weeks with 60% isocaloric fructose. 
Time 0 was considered as control (baseline) for each dog, providing the best reference for 
the progression of the metabolic syndrome. 

2.2. Hemodynamics 

Hemodynamic measurements were recorded in the conscious state using implanted 
instruments as previously described [24-27]. Coronary blood flow (CBF) and mean 
coronary flow (MCF) were determined with a Doppler flow transducer placed on the left 
circumflex coronary artery. The left ventricle was monitored for systolic (SBP) and 
diastolic blood pressure (DBP), mean arterial pressure (MAP) and end-diastolic (LVEP) 
pressure. The rate of the left ventricle pressure rise (dP/dt) was obtained with an 
operational amplifier (National Semiconductor LM 324). 

2.3. Echocardiography 

Echocardiography was performed with an Acusan 356 Sequoia in both a longitudinal 
and M-mode view [28] to determine the heart rate (HR), stroke volume (SV) and left 
ventricle diameter in diastole (LVDD) and systole (LVDS). Cardiac output (CO = HR x 
SV), fractional shortening (FS = (1-LVDS/LVDD)x100%), coronary vascular resistance 
(CVR = MAP/CBF) and total peripheral resistance (TPR = MAP/CO) were considered 
instead of the independent variables used to calculate them. 

2.4. Blood chemistry 

Blood samples were taken from an implanted catheter. Plasma insulin (INS), 
angiotensin II (ANG), uric acid (UA), homocysteine (HC), and high (HDL) and low 
(LDL) density lipids were measured as previously described [29]. 

2.5. Glucose tolerance test and insulin resistance 

In addition to dosing blood levels of glucose (GLU) and INS, we performed an 
intravenous glucose tolerance test in the morning, after 12 hours of fasting at baseline, 4 
and 7 weeks on fructose. Insulin resistance IR was defined as the insulin level [JNS] 
multiplied by the time necessary to reduce the glucose concentration at half, denoted as 
Tir. In normal conditions, increase of insulin level reduces T7,. 


IR =[INSJT,,, (1) 
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2.6. Data transformation 

Pathophysiological characteristics are of diverse nature, expressed in different units 
and bear unequal importance for the system, making difficult their integration into 
functional models. For instance, has 20% increase of SBP equivalent effects on the 
cardiovascular system as 20% increase of INS level, 20% increase of LDL or 20% 
increase of TPR? 

In order to make the alterations comparable, we transformed the real feature i values 


for dog 6 after w (= 0,1,2,3,4,5,6,7) fructose-feeding weeks, ow into numbers. First, we 
referred each feature a of each dog 6 at each time point w (weeks fructose feeding) to its 


baseline (self- reference to normalize the biological variability). Then we determined the 
mean and standard deviation for all dogs at the same time point. 


po = ia (0) — ( pie sg -4 
; q?) ps all 5 ia 
1 
g (2) 
ES : 
(w) ( cst 
Ss, = (a dips 
(5,0) . 
stdev(B, re if w>0 


To each feature i at each time point w we assigned the standard numerical score: 


ae w) 


jas or SoS a (3) 


In Gans foniation (3) we have implicitly assumed that the ratio-to-physiological 
feature values are from (or could be approximated with) a normal or a lognormal 
distribution [30] of possibilities. The lognormal distribution is a continuous probability 
distribution of a random variable whose logarithm is normally distributed. Almost all 
pathophysiological features spanning the entire axis of real numbers are normally 
distributed and almost all non-negative features are lognormally distributed. All features 
used in this report are lognormally distributed. The standard scoring (3) makes all 
lognormally distributed features comparable. Thus, a 20% increase of ¢ (but not of a) 
means the same deviation from normal, regardless the nature of that feature. 

2.7. Significant regulation 

A feature was considered as significantly regulated after w weeks of fructose feeding 
with respect to control if the as criterion (4) is satisfied: 


pooled coefficient of variation (4) 


Pra| test {a}, {a"},2,3) |< 0.05 
ENE Se ale 


heteroscedastic t-test of means equality 
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Criterion (4) requires that the absolute ¢ score exceeds the cut-off ¢ determined for 
that feature in that condition, AND the p-val of the heteroscedastic t-test for the equality 
of distributions’ means [31] is less than 0.05. The cut-off incorporates the technical noise 
and biological variability in both compared conditions (0 and w weeks of fructose 
feeding). 


2.8. Correlation of the pathophysiological characteristics 

One may consider the blood and the left ventricle of each of the six dogs, at each 
time point, as of one generic dog tuned to adapt six slightly different (not significantly 
regulating) local environments. Therefore, we have used Pearson product-moment 
correlation coefficient to determine the degree by which dog—to-dog variability in one 
feature is correlated to variability in each other feature in the same condition. Pearson 
correlation indicates whether the two features are linearly related. Also, we have 
determined the correlation between the profiles of two features for the same dog during 
the fructose feeding and averaged the results for the six dogs. In both cases, close to unit 
positive/negative values indicate synergistic/antagonistic expression, while close to zero 
values indicate independence of the considered features. The cut-off value of the 
coefficient of correlation between two features to be considered as synergistically/ 
antagonistically/independently expressed was determined for each feature pair by taking 
into account the effect of the inherent technical noise on those features [32]. However, 
only the correlation of dog-to-dog variation of two features in the same condition (# 
weeks of fructose feeding) indicates that such features may regulate each-other. 
Correlation between changes of two features during the seven-week experiment indicates 
only that a common factor (here fructose feeding) regulates both. 


2.9. The pre-Hilbert space of states and the “Patholog” of a state 

For quantitative evaluation of the fructose feeding effects, the model variables z are 
expressed to take positive/negative values for statistically significant up-/down- 
regulation, O-value indicating no regulation: 


gs 
(le -E*)) if a" > on 1 Poa < 0.05 


aK 
0 if else 


1 

Using the z-scores one can build the pre-Hilbert space (norm derived from a scalar 
product) of states having as many dimensions (NV) as independent features are considered 
from the perspective of Pearson correlation coefficient. In this space, each point 
represents a possible state and each curve a possible evolution. As constructed, the central 
point (all coordinates are 0) is the control state (here before starting the fructose diet). 
Given the independence and the way we normalize the features the space is orthonormal 
and the multi-dimensional Pythagoras’ theorem can be applied. Therefore, we define the 
“patholog” (sic!) Y of a state with respect to a subset of features (e.g. blood hormones) 
as the Euclidian distance to the control state [33, 34]. 


BOD) = f(z) (6) 


icD 


(5) 
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As defined, the “patholog” incorporates alterations of all quantified features in the 
selected subset. Since the feature subsets could be of different sizes (e.g.: 5 
hemodynamics vs 2 lipids) we defined and used the “standard patholog” IT 


= (2° y) (7) 


to compare the alterations of subsets of various sizes. 


2.10. Data processing 
OriginPro 8.5, GraphPad Prism 9 and Excel (Office 2020) were used for data analysis 


and graphs. 
3. Results and discussions 


3.1. Insulin resistance 
Figure 1 shows that insulin level and half-life in the glucose tolerance test are 
inversely correlated at the base line, justifying our IR index. Thus, normally fed dogs with 
higher insulin level metabolize faster (smaller half-life) the injected glucose than dogs 
with lower insulin levels. 


25 1.8 Fig. 1. Insulin level [INS] and half-life T12 in 

glucose tolerance test of normally fed dogs. 
Z 20 r The dogs were ordered with respect to 
=| S$ decreasing [INS]. Experimental data 
= " 3 were fitted with sigmoidal functions only 
S40 = to show that a decrease in the insulin 
‘2 level is accompanied by an increase of 
ri aa the glucose half-life (no relationship 

1 2 3 #46 5&5 6 between dog number and [INS] or T12). 


DOG# 


As presented in Fig. 2A, insulin level of high fructose-fed dogs increases 
exponentially (e.g. 3.4x after 4 weeks and 7.2x after 7 weeks) without exhibiting 
saturation during this period. The exponential growth constants were obtained using 
GraphPad Prism: 


INS(t)=INS(O)e“ => INS(w) = (10.47 £1.82)eO?P 0” 


However, instead of the expected faster glucose decay we found a substantially 
slower one, taking longer to reduce glucose concentration to half (i.e. larger half-life 
Tin, Fig. 2B). Experimental data were fitted with one phase decay law, with the decay 
rate K and [glu], and [,glu]_.. the initial (1 min after injection) and residual (at the end 
of the test) glucose levels: 


[slu], =(L2tu], [elu], Je" +[ elu], . half life T,, == 
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We found that /R increased faster than the insulin level (9.8x and 30.5x for IR vs 
3.4x and 7.2x for [INS]) . These results confirm the fructose-feeding induced insulin 
resistance syndrome [19, 35-39]. Using the above-defined IR measure, the insulin 
resistance increased by 9.75x after 4 weeks and by 30.45x after 7 weeks of fructose 
feeding (Fig. 2C). 


k 0 weeks (T4)2= 1.476 min) 609 


4 weeks (T4)2 = 4.760 min) 


“ ‘ a 2504 7 weeks (14/2 = 7.045 min) = 600 
E75 é a = 
2 d E = 
cE ar = 5 
= 50 Fl g 700 E 400 
3 72: 8 5 
e on Fy Fe 1504 3 195 
eT & 200 
: 1004 
20 
; — : : 6 
O42 6 @ 6 6 F 0 20 40 60 80 0 4 7 
weeks of high fructose diet minutes after administration weeks of high fructose diet 


Fig.2. Insulin level, glucose tolerance test and insulin resistance. A. Increase of plasma insulin 
level during high fructose feeding. B. Glucose test in 6 adult male dogs after 0, 4 and 7 weeks of 
fructose feeding. C. Increase of insulin resistance (IR). 


3.2. Alteration of hemodynamics, echocardiography, blood chemistry and lipids 

Figure 3 presents the evolution of several characteristics during high fructose diet. 
The absolute values were C-transformed to make comparable the severity of their 
alteration as illustrated in Table 1 for SBP and INS. 


al 


Baseline 


a. 

125 

120] 125|_100 
135) 
140 
50 


120) 140] 145] 135 
| weeks | 140| 160] 135] 140] 150] 140] 144.17] 9.17/ 
150.83] 8.01 
150 165 


Baseline E d 12.06 
Week 1 : 11.53 


41.05 
47.13|20 
64.80 


Table 1. Systolic Blood Pressure (SBP) and Insulin level (INS) real data and their 
transformation. Red background in column z indicates significant up-regulation. 


Note the significant increase of insulin (INS) and angiotensin II (ANG ID) after 3 
weeks and of “bad cholesterol” (LDL) and systolic blood pressure (SBS) after 4 weeks 
of high fructose feeding. Decrease of “good cholesterol” (HDL) and increase of 
homocystein (HC) and uric acid (UA) became significant after 5 weeks and increase of 
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diastolic blood pressure (DBP) after 6 weeks. The alterations continue to develop and 
the left ventricle end-diastolic (LVEP) pressure, coronary vascular (CVR) and total 
peripheral (TPR) changed significantly after 7 weeks. The heart rate, cardiac output, 
rate of left ventricle pressure rise, left ventricle diameter in systole/ diastole and the 
mean coronary flow were not significantly altered even after 7 weeks fructose. 
Interestingly, alteration of low-to high density lipids ratio (LDL/HDL) became 
significant after 2 weeks although the LDL changed significantly after 4 and HDL after 
5. 


A. HEMODYNAMICS B. ECHOCARDIOGRAPHY 


10 
cH) @ 
3 3 
m7 re 
we én 5 

0 a 
0 2 4 6weeks 0 2 4 6weeks 
C. CHEMISTRY D. LIPIDS 


— LDL 
—- HDL 
-e LDL/HDL 


c-score 


0 2 4 6weeks 0 2 4 6weeks 


Fig. 3. Evolution of hemodynamics, echocardiographic characteristics and levels of 
hormones and lipids during the 7 weeks high fructose feeding of 6 adult male dogs. 
Numbers after feature acronyms are the average + sd of the absolute control values (i.e. before 
starting the fructose diet). Data were normalized and transformed to numbers. Larger symbols 
indicate (p < 0.05) significant regulation with respect to control. Data were plotted on the same 
scale to compare the alterations of various features. 
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From the C-score perspective, increase of INS (¢ = 13.71), SBP (8.42) and 

LDL/HDL increase (7.53) had the highest impact on the dog health. However, none of 
these reported altered features exhibited any saturation and more alteration is expected 
if the fructose-feeding continued beyond the 7 weeks. 
We found that more features became altered when the duration of fructose feeding 
increases. Insulin (INS) had the largest alteration (C = 13.71) after 7 weeks fructose and 
MCF, dP/dt, CO and HR were not significantly altered at any time during this period. 
The difference between the levels of low (LDL) and high density (HDL) lipids was also 
plotted owing to its own clinical significance beyond the levels of the two components. 


3.3. Correlation of pathophysiological features 
No significant correlation was found between variations of MCF, SBP, DBP, LVEP, 
dP/dt, HR, CO, FS, CVR, TPR, INS, UA, HC, ANG, LDL and HDL in any condition. 
Therefore, all of these variables were considered as independent and the pre-Hilbert space 
of states was constructed with them. However, SBP, DBP, INS, UA, HC and LDL/HDL 
exhibited correlated increase in their absolute values during fructose-feeding (Fig. 4). 
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Fig.4. Correlated increase of hemodynamics, chemistry and low-to-high lipid density ratio 
during fructose feeding. p = pair-wise Pearson correlation coefficient between the average (for 
the six dogs) increase of the feature represented on x-axis and those of the features represented 
on the left and right axes. 0, 1, 2 , ... close to symbols indicate the number of weeks of fructose 
feeding. Note the significant correlated increase of these features. 


3.4. The pre-Hilbert space of pathological states 

Figure 5A presents a 3D section (systolic blood pressure, difference between 
the plasma levels low and high density lipids and insulin level) from the space of 
states. Figure 5B illustrates the 7-week dynamics of the patholog and Figure 5C 
presents the 7-week dynamics of the standard patholog. 
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Fig. 5. A. 3D section through the pre-Hilbert of states of a group of 6 adult male dogs 
subjected for 0, 1, ..., 7 weeks to 60% isocaloric fructose diet (note the projections in the 2D 
planes). Patholog for each time point with respect to these variables is the distance from the 
origin. B. Evolution of the patholog of all and subsets of features. Numbers indicate how 
many independent features were considered in each set. The plot started after the second week of 
fructose feeding since no pathologic alteration was recorded for one week. Note that although the 
Mets was installed after 4 weeks fructose feeding, the patholog more than doubled after other 3 
weeks fructose, without any indication of saturation. C. Evolution of the standard patholog for 
all and subset of features. Note that the largest alteration is associated with changes in 
hormones’ levels. 


Conclusions 


The dog was the first and it is still the most used large animal model for the human 
cardiovascular diseases (>54,000 published studies on canine models). However, dog’s 
heart is highly aerobic and very exercise tolerant, and metabolism evolution has favored 
protein diet as opposed to the humans which are omnivores. Although the Ossabaw pig 
model may be a better representation of genetically disposed human MetS [35-37], we 
have preferred the dog because it can be trained, allowing studies on previously 
instrumented, conscious animals with no anesthetic effects on the cardiovascular 
physiology [24-29]. 

We have observed that decrease of the insulin level is associated with increase of the 
half-life of the administrated glucose in the glucose tolerance test (Fig. 1). This expected 
observation led to a very simple and intuitive measure for the insulin resistance by 
linking the blood level of insulin to the decay of glucose level in the glucose tolerance 
test. Possible (not mutually exclusive) causes of the much faster increase of /R than of 
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[INS] include: i) saturation of insulin receptors, ii) reduction of insulin receptors or iii) 
alteration of insulin signaling pathway. An RNA-sequencing experiment is underway to 
check the last two potential causes at the transcriptomic level of the left ventricle. 
However, [INS] was not measured during the 75-min glucose tolerance test. Therefore, 
experimental data cannot be used in mechanistic models of insulin sensitivity in dogs 
subjected to intravenous glucose tolerance tests [38] or other “minimal models” [39] 
seeking fitting of experimental curves. 

Our mathematical approach with the “patholog” as a global quantifier offers the 
simplest, yet the most intuitive way to quantify and compare the effects of altering 
individual and combined features of wide diversity. Thus, if a feature is of higher 
importance, then the homeostatic mechanisms will keep its value (expression level for 
genes) within a narrow interval. Therefore, the c-transformation made comparable the 
effects of altered characteristics regardless of their different nature and measuring scale. 
For instance, 131.5 mmHg for SBP (i.e. an increase with 7% from the control 123 + 4 
mmHg in Fig. 2) has the same ¢ = 2 score as 29.0 uU/mL for INS (+141%), 195.9 mg/dL 
for HDL (+44%) or 2.0 mg/dL for UA (+70%). Moreover, the score ¢ is independent on 
the scale used for the absolute values: e.g. c = 2 for both 37.8°C (an increase with 2.2% 
from the normal 37.0 + 0.4°C) and 100.0°F (+ 1.0% from the normal or 98.6 + 0.7°F) in 
the case of orally taken body temperature. 

With the “patholog” values determined before and after a treatment (©) one can 
compute the efficacy F of that treatment as the percent reduction of the “patholog”: 


YF after 
E(®) =| 1-——— _|x100% (8) 
before © 
Such definition offers an objective score of the global amelioration of the disease- 
induced changes of physiological features. From this perspective, the possible outcomes 
are: 


1) E = 100% — oer = O (cured) 

2) 0< E < 100% —0 ¥< Paper < Yrefore (ameliorated) 
3) E=0 — Woper = Prefore (null) 

4) E<0< YWoer > Poefore (negative) 


Outcome (2) does not mean necessarily that the “after” state is identical to the 
“before” one. Instead, it means that the representative points are on the same 0-centered 
hypersphere (hence equal distance to the center, i.e. physiological state). 

In summary, our approach is the simplest (yet strong enough) one can imagine to 
globally quantify the evolution of a disease and result of a treatment. 


Notations 

CVR = coronary vascular resistance 
dp/dt = rate of LV pressure rise 

FS = fractional shortening 
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HDL = high density lipids (“good” cholesterol) 
HFD = high fructose diet 

IR = insulin resistance 

LAD = left anterior descending 

LDL = low density lipids (“bad” cholesterol) 
LV = left ventricle 

LVDD = left ventricle diameter in diastole 
LVDS = left ventricle diameter in systole 
MAP = mean arterial pressure 

MCF = mean coronary flow 

MetS = Metabolic syndrome 

RCA = right coronary artery 
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